8Ш 


УДК 519.651.5 


В.И. Шмойлов 
Южный научный центр РАН, г. Таганрог, Россия 
туз (@иге.ги 


О некоторых применениях 
парадоксального способа суммирования 
непрерывных дробей 


Рассматривается иное, нежели традиционное, определение сходимости непрерывных дробей. Новый 
метод суммирования используется при определении значений расходящихся в классическом смысле 
непрерывных дробей и рядов, которые нередко возникают при математическом моделировании тех 
или иных практически важных задач. 


Академик Н.Н. Моисеев в своей книге «Математика ставит эксперимент» [1] 
писал: «Нередко наше отставание в вычислительной технике нам удаётся компен- 
сировать более эффективными алгоритмами, применяемыми при решении задач». 
Вниманию читателей предлагается парадоксальный способ определения значений 
расходящихся непрерывных дробей, который может быть эффективно использован 
разработчиками прикладного программного обеспечения современных отечествен- 
ных вычислительных систем. 

Одним из факторов широкого использования непрерывных дробей в вычисли- 
тельной математике является то, что непрерывные дроби в большинстве случаев 
дают гораздо более общие представления трансцендентных функций, чем классические 
представления степенными рядами. Непрерывные дроби, зачастую, могут быть с 
большим эффектом использованы для ускорения сходимости медленно сходящихся 
рядов. Более того, преобразуя расходящиеся ряды в соответствующие непрерывные 
дроби, нередко можно просуммировать, то есть найти значения расходящихся рядов. 
Известно, что непрерывные дроби тесно связаны с аппроксимациями Паде, которые, 
как отмечается в [2], стали главным вычислительным средством в задачах статической 
механики и физики твердого тела. Поэтому существенные результаты, полученные в 
теории непрерывных дробей, в частности, в вопросах сходимости, могут быть исполь- 
зованы и в аппроксимациях Паде. 

Бесконечной цепной дробью, или непрерывной дробью, называют выражение вида 


а 
+ - ь 
Ь + 2 
+; а, 
1+ 
Ь +: 
где а; иБ,, {=1,2,... - в общем случае независимые переменные. 
Часто непрерывную дробь записывают в компактном виде в форме Гершеля: 
ы+ В ы 


В +Ь, Е ВЕ . 

Непрерывная дробь называется сходящейся, если последовательность ее 
подходящих дробей имеет конечный предел. Этот предел называется значением 
непрерывной дроби. Непрерывная дробь называется расходящейся, если последова- 


«Штучний 1нтелект» 42008 721 


Шмойлов В.И. 
8Ш 


тельность ее подходящих дробей предела не имеет. Имеется большое количество 
признаков сходимости, при помощи которых можно сказать, существует предел 
последовательности подходящих дробей или нет. Наиболее широкое применение, 
пожалуй, получил достаточный признак Ворпицкого [2]. По признаку Ворпицкого 
непрерывная дробь 

а @> а, 


1+14...+1+” 
<1/4, п=2,3.... 

В статье будет рассмотрено несколько задач из разных разделов вычисли- 
тельной математики, решенных при помощи так называемого ,/ф-алгоритма, — 
нового метода суммирования расходящихся непрерывных дробей. 

В [3] предложено иное, нежели традиционное, толкование сходимости не- 
прерывных дробей. Для установления значений непрерывных дробей будем 
использовать „/ф -алгоритм: 

Непрерывная дробь сходится и имеет своим значением в общем случае комп- 
лексное число ==ке% ‚, если существуют пределы 


сходится, если |а, 


5 


Вт ШВ О,|=%, (1) 
. К 

пит = ©, Ь (2) 
5—0 5 


где Р; /О; — значения 1-й подходящей дроби из совокупности, включающей 5 подхо- 
дящих дробей, К; — число отрицательных подходящих дробей из 5 подходящих 
дробей. 

Этот алгоритм применим как к обыкновенным непрерывным дробям, так и к 
непрерывным дробям с графами иных типов, например, к непрерывным дробям 
Хессенберга и ветвящимся непрерывным дробям. 

В случае непрерывных дробей, сходящихся в классическом смысле, аргумент 
ф, примет значения 0 илил. Если ф, = 0, то значение сходящейся непрерывной 
дроби будет совпадать со значением модуля их, : 

д = пе’ = п . 

Если ф, = л, то значение сходящейся непрерывной дроби будет отрицательное 

число: 


в=ке” =—к. 
Предложенный г/ф-алгоритм даёт возможность устанавливать значения расхо- 


дящихся непрерывных дробей. Проиллюстрируем эффективность предложенного 
способа суммирования расходящихся в классическом смысле непрерывных дробей 
решением ряда задач вычислительной математики. 

Известно, что непрерывные дроби целесообразнее использовать для аппрок- 
симации функций, нежели степенные ряды, так как непрерывные дроби зачастую 


сходятся в более широкой области. Например, ряд Меркатора 
2 3 4 
х 
№ +х) =х-—+ 
РЕ 59 —.3 
представляет логарифмическую функцию в единичном круге, в то время как 


непрерывная дробь Лагранжа 


ххх 2х 2х пх пх 
Ш(+х) = 
Т+2+3+2 +5 +...+ 2 +21 41+... 
сходится к функции ш([+х) на всей плоскости комплексного переменного, за 


исключением выреза от —1 до — [2]. 
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При отрицательных значениях аргумента логарифмическая функция имеет комп- 
лексное значение, и естественно, что непрерывные дроби, получающиеся из непрерывной 
дроби Лагранжа при х<-1 будут расходящимися в классическом смысле. 

На рис. 1 (а, 6, в, г, д, е) показано распределение подходящих непрерывной 


дроби Лагранжа (2.1) при х = -10, —100, —1000 на начальных участках 
(п=1-100 и п=1-500). 


(9) = 3.8337е' 78094 


(9) = 3.8337е' 96094 


№1(-99) = 5.5663е 


600 = 5.5663" 79 


10.4268 


11(-999) = 7.5876е 


Рисунок 1 — Распределение значений подходящих непрерывной дроби 
логарифмической функции 
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В табл. 1 показаны результаты суммирования расходящейся непрерывной дроби 
1-2-3-2-5 2-2п+1 
При отрицательном аргументе логарифмическая функция имеет комплексное 
значение 


№ (2) = 


№ (—2) = 3,2171505117е/' 5363945 , 
которое, естественно, не может приближаться непрерывной дробью с вещественными 
элементами и, тем не менее, суммирование при ПОМОЩИ "/ ф -алгоритма позволяет 


установить значение дроби (3). 


Таблица 1 


16 
32 
64 
128 
256 
512 
1024 
2048 
4096 
8192 
16384 
32768 
65536 
131072 
262144 
524288 
1048576 
2097152 
4194304 
8388608 


Значение 
подходящей 


—3.0000000 
6.0000000 
—3.0000000 
—97.5000000 
1.4880473 
3.1985122 
62.8693924 
0.9165216 
1.7095765 
3.9037050 
—15.4772571 
2.6358581 
11.1007665 
—0.6961262 
—1.7591587 
—6.4347291 
5.5879135 
—3.9038315 
16.0431708 
—0.0551483 
—0.2709104 
—0.7308612 
—1.8537413 
—7.2124648 


Модуль 
комплексного 
числа, Г, 


3.0000000000 
4.2426406871 
3.0000000000 
4.9614481602 
3.5474336503 
3.6050160485 
3.3885474566 
3.1810462758 
3.2148854739 
3.2112688498 
3.2219262392 
3.2194825453 
3.2127253440 
3.2169015620 
3.2167104407 
3.2170964982 
3.2171496506 
3.2171884212 
3.2171480639 
3.2171287791 
3.2171427009 
3.2171496552 
3.2171502478 
3.2171495794 


Погрешность, 


=, = —", 


0.2171505117 
1.0254901754 
0.2171505117 
1.7442976485 
0.3302831386 
0.3878655367 
0.1713969449 
0.0361042359 
0.0022650377 
0.0058816618 
0.0047757275 
0.0023320336 
0.0044251676 
0.0002489496 
0.0004400709 
0.0000540134 
0.0000008610 
0.0000379094 
0.0000024477 
0.0000217325 
0.0000078107 
0.0000008564 
0.0000002638 
0.0000009323 


3.1415926535 
1.5707963267 
1.5707963267 
1.5707963267 
1.3744467859 
1.3744467859 
1.3744467859 
1.3499030933 
1.3621749396 
1.3499030933 
1.3560390164 
1.3529710549 
1.3529710549 
1.3533545501 
1.3533545501 
1.3536421715 
1.3536421715 
1.3536182030 
1.3535942346 
1.3536421715 
1.3536361793 
1.3536361793 
1.3536391754 
1.3536399244 


Определение значения расходящейся непрерывной дроби 


Погрешность, 
бо т. | = Ф, 


1.7879528081 
0.2171564813 
0.2171564813 
0.2171564813 
0.0208069405 
0.0208069405 
0.0208069405 
0.0037367521 
0.0085350941 
0.0037367521 
0.0023991710 
0.0006687905 
0.0006687905 
0.0002852953 
0.0002852953 
0.0000023260 
0.0000023260 
0.0000216423 
0.0000456108 
0.0000023260 
0.0000036660 
0.0000036660 
0.0000006699 
0.0000000790 


33366 3п Зп 
23255 2-2п+1 
7, =3.2171505117... Фо =1.3536398454.. 


В первой колонке таблицы даны номера и подходящих дробей разложения (3). 


№ (2) = 


Номера подходящих дробей составляют степень 2: и=2',{=1-+23. Значения подходя- 
щих дробей с этими номерами приведены в соседней колонке 2. Как и следовало 
ожидать, значения подходящих дробей {Р‚ /О, | с ростом и не стремятся к какому-либо 
пределу. Для чисел же, расположенных в колонке 3, напротив, стремление к пределу 
можно без труда обнаружить — значения асимптотически приближаются к величине 
3.2171505117..., то есть к модулю комплексного числа ш(-2). Даже беглого взгляда 
на колонки 6 и 7 достаточно, чтобы убедиться, что с ростом количества подходящих 
дробей разложения (3) все более точно устанавливается значение аргумента иско- 
мого комплексного числа. 
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Найдены представления элементарных и некоторых специальных функций в 
виде непрерывных дробей Хессенберга [4]. Например, 


х/И х/2 х/3! х/4Й х/5! 
РЕ ХИ х/21 х/ЗЕ х/4 
оо фт ИИ 2/21 д 
оо и х/21 
ооо и... 
1 а. | . . НЙ. (4) 

и (1+ х/И х/2! х/3З! х/4 ... 
Го х/Ш х/21 х/З 

оо фр х/21 
оо фри 


Непрерывная дробь (4) определяет логарифмическую функцию на всей плоскости 
комплексного переменного с разрезом от 0 до —1. Комплексное значение логарифми- 
ческой функции на разрезе определяется из (4) суммированием по формулам (1) и (2). 
В качестве подходящих дробей (4) будем брать последовательность отношений опре- 
делителей: 


х/И х/21 х/З 
х/! = То ХИ жи 

рей. Е Е ЭН о 
1 х/1! ь х/И х/2|” 


Е о х/Н 
В табл. 2 приведены результаты вычисления комплексного числа 1/(-2) при 
помощи функциональной непрерывной дроби Хессенберга (4) и алгоритма суммиро- 
вания расходящихся непрерывных дробей, то есть формул (1) и (2). 
Определение значения непрерывной дроби (2.3) при х=-\З. 
х=-И3, п =0.3108340739.., ф,=-1.3536398454.. 


Таблица 2 


Модуль Аргумент 
Значение У. р рум о 
подходящей |комплексного ши | комплексного =, В о пт 
и ый и, числа, ==” на 


Нее Нова 

0.82455439 0.51372032 

0.50977666 0.19894259 

0.40513373 0.09429965 

0.35654051 0.04570644 

0.33285231 0.02201823 

0.32176153 0.01092745 

0.31569144 | 0.00485737 

0.31342156 0.00258748 

. 0.31217034 0.00133626 
0.30407891| 0.31149855 0.00066447 
1.28357034| 0.31116989 0.00033582 
—0.09454281| 0.31099243 0.00015836 
0.31091665 0.00008257 

0.31087600 0.00004193 

0.31085447 0.00002039 

0.31084454 | 0.00001047 

0.31083902 0.00000495 

0.31083657 0.00000249 

. 0.31083533 0.00000126 
11010048 . 0.31083471 0.00000064 


НЕ 41371669 Ноя — 600768 
—1.34639685 0.00724299 
—1.42119667 0.06755683 
—1.38379676 0.03015691 
—1.36509680 0.01145696 
—1.35574682 0.00210698 
—1.35574682 0.00210698 
—1.35574682 0.00210698 
—1.35457808 0.00093823 
—1.35399370 0.00035386 
—1.35370152 0.00006167 
-1.35370152 0.00006167 
—1.35370152 0.00006167 
—1.35366499 0.00002515 
—1.35364673 0.00000689 
—1.35364673 0.00000689 
—1.35364217 0.00000232 
—1.35363988 0.00000004 
—1.35363988 0.00000004 
—1.35363988 0.00000004 
—1.35363988 0.00000004 


ЕЕ 
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Приведём ещё пример нахождения значений расходящейся в классическом 
смысле непрерывной дроби при помощи описанного выше алгоритма суммирования. 

Известна непрерывная дробь Лагранжа для степенной функции [2]: 
ух (1-у)х ЦИ+и)х (т-у)х  (т-+у)х (5) 
1 2 жа 2+ 2т+1+... 

Непрерывная дробь (5) сходится на всей плоскости комплексного переменного, 
разрезанной по вещественной оси от-1 до —©. При х<-1 подкоренное выражение 
функции у=(1+х)” становится отрицательным, следовательно, значение функции (1+х)” 
при х<-1 становится комплексной величиной, которая, очевидно, не может прибли- 
жаться вещественной последовательностью подходящих дробей разложения (5), 
элементы которого при х<-1 остаются действительными величинами. 

В табл. 3 приведены результаты вычисления значения расходящейся непрерыв- 


(1+х)” =1+ 


ной дроби для \/-9 при помощи „/ф -алгоритма. 


Определение значения расходящейся непрерывной дроби 
9-10 240 4410 540 (3п-10:10 (3и+1-10 
3- 2 9 5 5 = ВО. 


’, =2.080083..., 
Таблица 3 
Значение 


подходящей 
дроби 
2 


Номер 

звена 

дроби 
1 


2.42857 14 
7.9230769 
0.0083609 
2.3793179 
3.4334933 
—6.6957404 
1.8056646 
2.7046939 
14.0355053 
0.4131559 
—0.5703198 
—33.5321135 
1.1197769 
1.0914866 
1.0349376 
0.9216825 
0.6914169 
0.1869493 
—1.4500205 
5.5741164 
—0.8763266 
15.7595707 
0.4758782 


дроби для \/-9. 


Фо =1.047197.... 


Модуль 
комплексного 
СИЕ ", 


ОЗ 0952 8095238 
1.892578548704 
0.94690843 1432 
1.514212477498 
1.931910874242 
2.285828469458 
2.171608471371 
2.095947975804 
2.077080542151 
2.063579787259 
2.085479528230 
2.082858835017 
2.079900863342 
2.080396700069 
2.080307362394 
2.079865087605 
2.080106476321 
2.079930141859 
2.080071498760 
2.080081877069 
2.080075907132 
2.080077041783 
2.080083260804 


Аргумент 
и шш | комплексного 
“== НЕ числа, я 


ТЕТЕ ЕЕ 
0.187505274347 
1.133175391618 
0.565871345553 
0.148172948809 
0.205744646406 
0.091524648319 
0.015864152753 
0.003003280900 
0.016504035792 
0.005395705178 
0.002775011965 
0.000182959709 
0.000312877017 
0.000223539343 
0.000218735446 
0.000022653270 
0.000153681192 
0.000012324291 
0.000001945982 
0.000007915919 
0.000006781268 
0.000000562247 


оо 
0.000000000000 
0.785398163397 
0.785398163397 
0.687223392972 
0.932660319034 
0.957204011640 
1.043106935762 
1.043106935762 
1.046174897338 
1.044640916550 
1.046174897338 
1.048092373322 
1.047133635330 
1.047229509129 
1.047229509129 
1.047205540679 
1.047193556454 
1.047199548567 
1.047190560398 
1.047195054483 
1.047199548567 
1.047197301525 


и" 
2 Ель Ф, | 


НТБ 55 196 Н- 


1.047197551196 
0.261799387799 
0.261799387799 
0.359974158223 
0.114537232162 
0.089993539555 
0.004090615434 
0.004090615434 
0.001022653858 
0.002556634646 
0.001022653858 
0.000894822126 
0.000063915866 
0.000031957933 
0.000031957933 
0.000007989483 
0.000003994741 
0.000001997370 
0.000006990797 
0.000002496713 
0.000001997370 
0.000000249671 


Рисунок 2 — Распределение значений подходящих непрерывной дроби для \/-9 
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Следует отметить, что при использовании г/ф-алгоритма необходимо иметь 


эффективные алгоритмы для построения так называемых соответствующих непрерыв- 
ных дробей и для вычисления длинных серий значений подходящих дробей. Исполь- 
зование для счёта классического рекуррентного алгоритма (ЕВ-алгоритм) приводит к 
быстрому переполнению разрядной сетки компьютера, а применение естественной 
процедуры вычисления непрерывной дроби «снизу-вверх» (ВВ-алгоритм) невозможно 
из-за недопустимо больших временных затрат при определении серий значений подхо- 
дящих дробей. В [4] были детально рассмотрены алгоритмы для вычисления длинных 
серий значений подходящих дробей. В [5] изложены различные методы построения 
соответствующих непрерывных дробей. 

Расходящиеся ряды нередко суммируются через соответствующие цепные 
дроби [6]. Для степенного ряда 


со +с1х+ вх КНС, №: (6) 
можно построить непрерывную дробь 
ях @х вх 9х @х 
©, + 1 2 3 21-1 2п р (7) 
о о в 


такую, что разложение п-й подходящей дроби этой цепной дроби будет совпадать с 


исходным рядом (6) вплоть до члена с„Х” включительно: 
Р,(х) 
0, (<) 


Такую непрерывную дробь называют соответствующей ряду или соответст- 
вующей непрерывной дробью. 


И 2 п п-1 
= С +ах+с,х +... Нс,Х НУ, их Я 


Используя коэффициенты с, степенного ряда (6), можно построить соответствую- 
щую непрерывную дробь (7) по формулам Хейлерманна — Стилтьеса [7]: 
0% = С,  =С,, 


@5„ = ФУ ? @„=- Фин У ; 
Ф, У, Ф, "Ули 
и С Ст С С Ст 
а С С .-. С та С 4 хх ба | (8) 
Сп Син Си Сп Син -**  С2и-2 
Фо =Ь Ч =Г. 
Найдем значение расходящегося ряда Эйлера: 
11122 
1— 142-341... = п” =0,5963473623231945... (9) 


1+1+1+1+1+...+1+1+... 


Ряд Эйлера (9) связан с интегральной показательной функцией ЕЦх). 
Действительно, известен быстро сходящийся ряд: 


—еЕ-1) = -4с-н ] = 0,5963473623231945... , 


= Ри. 
2.2! 3.3 


где е- неперово число, равное 2.718281..., 
С — постоянная Эйлера, имеющая значение 0.577215... 
Определив по расходящимся рядам соответствующие непрерывные дроби, 
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можем найти значения других расходящихся рядов. Например, 


1+1—1+2—5+14—42+132—.=Е Е 1 -1+ 5 | 618033... (10) 
1+1+...+1+... 2 
фе. 3 п 

Е ое = 0,655679 ... (11) 
ТЕ. 


Расходящийся ряд (10) имеет своим значением «отношение золотого сечения», 
а расходящийся ряд (11) связан с неполной гамма-функцией 


ро - [г а 66556795404. 
ЕЕ. 1295 2 и 


Метод суммирования расходящихся непрерывных дробей может быть эффек- 
тивно использован при решении такой практически важной проблемы, как решение 
бесконечных нерегулярных систем линейных алгебраических уравнений. Известно, 
что решение трёхдиагональной системы линейных алгебраических уравнений можно 
представить в виде непрерывных дробей, частными числителями и знаменателями 
которых были бы некоторые выражения из элементов исходной матрицы. Этот способ 
записи решений трёхдиагональной системы эквивалентен алгоритму «прогонки». 
Непрерывные дроби, представляющие решения системы, будут бесконечными, если 
бесконечна система линейных алгебраических уравнений. Но бесконечные непре- 
рывные дроби могут быть сходящимися и расходящимися. Задачи математической 
физики зачастую описываются бесконечными системами линейных алгебраических 
уравнений. Нередко при такой аппроксимации получают так называемые нерегулярные 
системы алгебраических уравнений. Под нерегулярными понимают системы алгеб- 
раических уравнений, решения которых не стремятся к пределу в классическом 
смысле с ростом размерности системы. 

При решении регулярных систем не возникает трудностей, так как вместо беско- 
нечных систем оперируют «усечёнными» системами. Непосредственно метод «усечения», 
очевидно, не срабатывает в случае нерегулярных бесконечных систем. При некоторых 
ограничениях решения нерегулярных бесконечных систем находятся при помощи 
рассмотренного выше метода суммирования расходящихся непрерывных дробей. Можно 
указать ещё ряд «проблемных» задач в вычислительной математике, эффективное ре- 
шение которых обеспечивается алгоритмом суммирования расходящихся непрерывных 
дробей. 
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Г.1. 5йтоуу 

Оп Зоте АррИсайоп$ 0Ё Рагадоха1 Мео4д о? Сопйпиед Егасйоп$ Зититайоп 

ш 15 агае \ме дезсге апофег, еп га@!юопа|, дейп1Аюп оЁ сопуегоепсе оЁ сопипие4 НасНопз. Ме\у 
тейо4 о# заттайолп 1$ изе4 Рог деегттайпе уа[аез оЁ уегоепе ш с1аз$1с зепзе сопипиеа Насйоп$ ап4 
зетез, умсЬ овеп арреагз Фигте Фе таетайса| тодейт» Фезе ог фозе ргасйсаПу ппроцапе ргоЫеттз. 
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